The  Astrophysical  Journal  Letters,  713:L136-L139,  2010  April  20 

©  2010.  The  American  Astronomical  Society.  All  rights  reserved.  Printed  in  the  U.S.A. 


doi:  10. 1088/2041-8205/7 13/2/Ll  36 


KEPLER-6b:  A  TRANSITING  HOT  JUPITER  ORBITING  A  METAL-RICH  STAR* 

Edward  W.  Dunham1,  William  J.  Borucki2,  David  G.  Koch2,  Natalie  M.  Batalha3,  Lars  A.  Buchhave4,5, 
Timothy  M.  Brown6,  Douglas  A.  Caldwell7,  William  D.  Cochran8,  Michael  Endl8,  Debra  Fischer9,  Gabor  Furesz4, 
Thomas  N.  Gautier  III10,  John  C.  Geary4,  Ronald  L.  Gilliland11,  Alan  Gould12,  Steve  B.  Howell13,  Jon  M.  Jenkins7, 
Hans  Kjeldsen14,  David  W.  Latham4,  Jack  J.  Lissauer2,  Geoffrey  W.  Marcy15,  S0ren  Meibom4,  David  G.  Monet16, 

Jason  F.  Rowe717,  and  Dimitar  D.  Sasselov4 
1  Lowell  Observatory,  Flagstaff,  AZ  86001,  USA 
2  NASA  Ames  Research  Center,  Moffett  Field,  CA  94035,  USA 
3  San  Jose  State  University,  San  Jose,  CA  95192,  USA 
4  Harvard-Smithsonian  Center  for  Astrophysics,  Cambridge,  MA  02138,  USA 
5  Niels  Bohr  Institute,  Copenhagen  University,  DK-2100  Copenhagen,  Denmark 
6  Las  Cumbres  Observatory  Global  Telescope,  Goleta,  CA  93117,  USA 
7  SETI  Institute,  Mountain  View,  CA  94043,  USA 

8  University  of  Texas,  Austin,  TX  78712,  USA 

9  Yale  University,  New  Haven,  CT  06510,  USA 

10  Jet  Propulsion  Laboratory /California  Institute  of  Technology,  Pasadena,  CA  91 109,  USA 
11  Space  Telescope  Science  Institute,  Baltimore,  MD  21218,  USA 
12  Lawrence  Hall  of  Science,  Berkeley,  CA  94720,  USA 
13  National  Optical  Astronomy  Observatory,  Tucson,  AZ  85719,  USA 
14  University  of  Aarhus,  Aarhus,  Denmark 
15  University  of  California,  Berkeley,  Berkeley,  CA  94720,  USA 
16  US  Naval  Observatory,  Flagstaff  Station,  Flagstaff,  AZ  86001,  USA 
Received  2009  November  15;  accepted  2010  January  8;  published  2010  March  30 


ABSTRACT 

We  announce  the  discovery  of  Kepler-6b,  a  transiting  hot  Jupiter  orbiting  a  star  with  unusually  high  metallicity, 
[Fe/H]  =  +0.34  ±  0.04.  The  planet’s  mass  is  about  2/3  that  of  Jupiter,  MP  =  0.67  Mj,  and  the  radius  is  30%  larger 
than  that  of  Jupiter,  RP  —  1.32  Rj,  resulting  in  a  density  of  pP  =  0.35  gem-3,  a  fairly  typical  value  for  such  a 
planet.  The  orbital  period  is  P  —  3.235  days.  The  host  star  is  both  more  massive  than  the  Sun,  M„  =  1 .21  M0,  and 
larger  than  the  Sun,  R*  =  1.39  R0. 

Key  words:  planetary  systems  -  stars:  individual  (Kepler-6,  KIC  10874614,  2MASS  19472094+4814238)  - 
techniques:  spectroscopic 


1.  INTRODUCTION 

The  Kepler  Mission  was  launched  on  2009  March  6,  to 
undertake  a  search  for  Earth-size  planets  orbiting  in  the  habitable 
zones  of  stars  similar  to  the  Sun.  Kepler  uses  the  transit 
photometry  approach  for  this  task  (Borucki  et  al.  2010a). 
Kepler’s  commissioning  process  went  very  well  and  the  system 
provides  data  of  exceptional  photometric  quality  (Koch  et  al. 
2010a).  Indeed,  the  final  commissioning  data,  9.7  days  of 
science-like  observations  of  50,000  stars  selected  from  the 
Kepler  Input  Catalog  (KIC;  Koch  et  al.  2010a),  were  of  such  high 
quality  that  they  are  now  commonly,  if  incorrectly,  referred  to 
as  the  “quarter  0”  data.  A  number  of  Kepler  Objects  of  Interest 
(KOIs)  simply  fell  out  of  this  data  set,  including  the  object 
discussed  here,  which  is  known  as  KOI-17.  Its  transit  signal  is 
huge  by  Kepler  standards,  and  excellent  light  curve  results  were 
obtained  quickly.  The  follow-up  spectroscopic  observations 
took  somewhat  longer  since  observations  were  not  planned  to 
be  undertaken  so  soon  after  commissioning.  As  a  result,  the  first 
data  intended  for  science,  33.5  days  of  observations  of  150,000 
KIC  stars,  became  available  in  the  meantime  and  were  folded 
in  with  the  original  test  data  for  the  analysis  presented  here. 
The  1.3  day  gap  between  these  two  sets  of  data  contained  one 


*  Based  in  part  on  observations  obtained  at  the  W.  M.  Keck  Observatory, 
which  is  operated  by  the  University  of  California  and  the  California  Institute  of 
Technology. 
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transit  of  Kepler-6b.  Haas  et  al.  (2010)  provide  additional  details 
regarding  the  Kepler  data  used  in  this  analysis. 

Because  the  pace  of  early  Kepler  extrasolar  planet  discoveries 
is  limited  by  the  radial  velocity  follow-up  observations,  the 
first  objects  to  be  announced  tend  to  be  those  most  easily 
confirmed  spectroscopically.  Kepler-6b  (=  KIC  10874614,  a  = 
19h47m20!94,  5  =  +48°14'23"8,  J2000,  KIC  r  =  13.253  mag) 
is  one  of  these  objects. 

2.  KEPLER  PHOTOMETRY 

The  Kepler  science  data  for  the  primary  transit  search  mission 
are  the  long  cadence  data  (Jenkins  et  al.  2010a).  These  consist  of 
sums  very  nearly  30  minutes  in  length  (with  rigorously  the  same 
integration  time)  of  each  pixel  in  the  aperture  containing  the 
target  star  in  question.  These  data  proceed  through  an  analysis 
pipeline  to  produce  corrected  pixel  data,  then  simple  unweighted 
aperture  photometry  sums  are  formed  to  produce  a  photometric 
time  series  for  each  object  (Jenkins  et  al.  2010b).  The  many 
thousands  of  photometric  time  series  are  then  processed  by  the 
transiting  planet  search  (TPS)  pipeline  element  (Jenkins  et  al. 
2010b).  The  candidate  transit  events  identified  by  TPS  were  then 
vetted  by  visual  inspection. 

The  light  curves  produced  by  the  photometry  pipeline  tend 
to  show  drifts  due  to  an  extremely  small,  slow  focus  change 
(Jenkins  et  al.  2010a),  and  there  are  also  sometimes  low- 
frequency  variations  in  the  stellar  signal  that  can  make  analysis 
of  the  transit  somewhat  problematic.  These  low-frequency 
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Time  HJD-2454900 
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Figure  1.  Detrended  light  curve  for  Kepler-6.  The  time  series  for  the  entire  data 
set  is  plotted  in  the  upper  panel.  The  lower  panel  shows  the  photometry  folded 
by  the  3.234723-day  period.  The  model  fit  to  the  primary  transit  is  overplotted 
in  red  (vertical  axis  on  the  left),  and  our  attempt  to  fit  a  corresponding  secondary 
eclipse  is  shown  in  green  with  the  expanded  and  offset  scale  on  the  right. 

effects  can  be  removed  by  modest  filters  that  have  only  an 
insignificant  effect  on  the  transit  signal  (Koch  et  al.  2010b). 
The  unfolded  and  folded  light  curves  for  Kepler-6b  produced  in 
this  manner  are  shown  in  Figure  l.18 

Figure  1  illustrates  two  of  the  ways  we  reject  false  positives 
based  on  the  Kepler  light  curves  alone.  The  lower  part  of  the 
figure  shows  the  planetary  transit  (referred  to  the  left  ordinate). 
The  points  marked  with  +  and  x  symbols  refer  to  alternating 
even  and  odd  transit  events.  This  aids  in  detecting  eclipsing 
binary  stars  with  nearly,  but  not  quite  equal  primary  and 
secondary  minima.  No  such  deviation  is  observed  for  Kepler-6. 
The  upper  curve  shows  the  time  of  secondary  minimum, 
assuming  zero  eccentricity.  The  vertical  scale,  on  the  right,  is 
magnified  and  measured  in  parts  per  million  to  help  search  for 
a  small  secondary  minimum  that  might  be  due  to  occultation  of 
an  M  dwarf  companion.  Again,  no  such  signature  is  visible. 

3.  FOLLOW-UP  OBSERVATIONS 

The  crucial  importance  of  ground-based  follow-up  obser¬ 
vations  was  recognized  long  ago  (Latham  2003)  and  a  well- 
established  plan  was  in  place  by  the  time  of  launch.  This  has 
evolved  somewhat  and  our  current  approach  is  described  by 
Gautier  et  al.  (2010).  The  key  features  are:  (1)  carry  out  high- 
resolution  reconnaissance  spectroscopy  to  determine  the  stellar 
properties  and  the  star’s  suitability  for  radial  velocity  work  and 
to  search  for  signs  of  stellar  companions,  (2)  search  for  faint 
blended  eclipsing  binaries  through  the  centroid  shift  approach 
(Batalha  et  al.  2010)  and  by  means  of  high-resolution  imag¬ 
ing,  and  (3)  detect  the  stellar  reflex  motion  with  precise  radial 
velocity  measurements.  In  the  future,  we  hope  to  improve  our 
constraints  on  stellar  properties  by  including  both  asteroseis- 
mology  results  for  at  least  the  brighter  host  stars  (Gilliland  et  al. 
2010a)  and  parallaxes  (Monet  et  al.  2010)  for  every  target.  These 
results  are  not  available  for  Kepler-6  at  this  time. 

18  The  time  series  photometry  and  radial  velocity  data,  including  bisector 
spans,  may  be  retrieved  from  the  MAST /HLSP  data  archive  at 
http  ://archive .  stsci .  edu/prepds/kepler Jtlsp . 


Table  1 

Relative  Radial  Velocity  Measurements  of  Kepler-6 


BJD 

Phase 

RV 

(ms-1) 

crv 

(ms-1) 

2454986.087289 

0.769 

+68.42 

4.48 

2454987.087348 

0.078 

-59.94 

4.59 

2454988.077695 

0.384 

-67.64 

4.75 

2454989.019649 

0.676 

+48.64 

4.64 

2454995.112273 

0.559 

+18.56 

4.79 

2455014.881155 

0.671 

+51.75 

4.54 

2455015.984052 

0.012 

-34.79 

4.65 

2455017.067098 

0.346 

-75.84 

4.94 

2455044.019956 

0.679 

+50.59 

5.42 

3.1.  Reconnaissance  Spectroscopy 

The  reconnaissance  spectra  for  Kepler-6  were  carried  out  very 
early  in  the  follow-up  program,  so  observations  were  obtained 
with  several  of  our  resources  for  comparison  purposes:  the  TRES 
spectrograph  on  the  Tillinghast  telescope  (by  L.  Buchhave  and 
D.  Latham),  the  McDonald  2.7  m  Coude  spectrograph  (by  M. 
Endl  and  W.  Cochran),  and  the  FIES  spectrograph  on  the  Nordic 
Optical  Telescope  (by  L.  Buchhave  and  G.  Furesz).  The  spectra 
indicated  that  Kepler-6  was  well  suited  for  precise  radial  velocity 
work  and  no  evidence  of  double  lines  or  a  stellar  companion  was 
seen.  The  best  data  on  stellar  properties  were  determined  later 
and  are  discussed  below. 

The  mean  radial  velocity  of  Kepler-6  was  determined  by 
C.  Chubak  with  two  High  Resolution  Echelle  Spectrometer 
(HIRES)  spectra  using  telluric  lines  to  set  their  zero  point. 
Observations  of  the  radial  velocity  standard  HD  182488  using 
the  same  approach  produced  a  mean  value  of  —21.50  ± 
0.18  kms-1,  in  good  agreement  with  the  adopted  value  of 
—21.508  kms-1.  Observations  with  FIES  and  TRES  using  the 
same  standard  star  to  determine  their  zero  points  were  in  good 
agreement  with  the  HIRES  velocity.  The  mean  radial  velocity 
for  Kepler-6  is  —49.14  ±  0.10  km  s-1. 

3.2.  High-resolution  Imaging  and  Centroid  Motion 

Kepler-6  has  a  companion  4"1  distant  and  3.8  mag  fainter 
as  seen  in  both  a  HIRES  guide  camera  image  and  a  PHARO 
J  band  AO  image  obtained  with  the  Palomar  Hale  telescope 
(Gautier  and  Ciardi).  A  WIYN  speckle  image  (Howell)  shows 
only  a  single  star  within  its  2"  square  box.  There  are  apparently 
two  additional  components  both  approximately  5.4  mag  fainter 
and  1  IT 5  distant  that  appear  in  the  KIC  as  KIC  IDs  10874615 
and  10874616.  These  two  objects  refer  to  a  single  star  and 
erroneously  appear  as  two  objects  in  the  KIC.  This  star  is 
very  close  to  the  edge  of  the  aperture  containing  Kepler-6.  The 
observed  shift  in  the  Kepler  intensity-weighted  mean  image 
centroid  during  transit,  (—0.1,  +0.3)  millipixels  in  the  (column, 
row)  directions,  is  consistent  with  the  expected  shift  if  the  faint 
companions  are  not  the  source  of  the  transit  signature.  If  the 
fainter,  more  distant  star  is  fully  in  the  aperture  the  expected 
shift  is  (0.0,  +0.3)  millipixels  and  if  it  is  completely  out  of 
the  aperture  the  expected  shift  is  (—0.1,  +0.4)  millipixels,  both 
being  consistent  with  the  observed  centroid  shift.  On  the  other 
hand,  if  the  closer,  brighter  background  star  is  the  transit  source 
the  expected  shift  is  30  times  larger,  ruling  out  this  possibility. 
The  fainter  object  is  too  faint  to  produce  the  observed  transit 
depth  even  if  it  disappears  entirely  during  the  transit,  and  if  it 
could  the  resulting  centroid  shift  would  be  nearly  a  factor  of  100 
larger  than  the  observed  shift.  The  typical  sensitivity  level  for 
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Table  2 

System  Parameters  for  Kepler-6 


DUNHAM  ET  AL. 


Vol.713 


Parameter 

Value 

Notes 

Transit  and  orbital  parameters 

Orbital  period  P  (d) 

3.234723  ±  0.000017 

A 

Midtransit  time  E  (HJD) 

2454954.48636  ±  0.00014 

A 

Scaled  semimajor  axis  «//?* 

7  OS"1-0'11 
'■u  - 0.06 

A 

Scaled  planet  radius  Rp/R* 

0-098292“ 

A 

Impact  parameter  b  =  a  cos  i  /  R+ 

0.398+„°o0023°9 

A 

Orbital  inclination  i 

86?  8  ±  0.3 

A 

Orbital  semi-amplitude  K(ms-1) 

80.9  ±  2.6 

A,  B 

Orbital  eccentricity  e 

0  (adopted) 

A,  B 

Center-of-mass  velocity  y  (ms-1) 

-18.3  ±  3.5 

A,  B 

Observed  stellar  parameters 

Effective  temperature  Teff  (K) 

5647  ±  44 

C 

Spectroscopic  gravity  log  g  (cgs) 

Fixed  at  transit  value 

C 

Metallicity  [Fe/H] 

+0.34  ±  0.04 

C 

Projected  rotation  v  sin  i  (kms-1) 

3.0  ±  1.0 

C 

Mean  radial  velocity  (kms-1) 

-49.14  ±  0.10 

B 

Derived  stellar  parameters 

Mass  M„  (Mg) 

1  90Q+0  044 

0.038 

C,D 

Radius  R*  ( Rq ) 

i  on  1+0.017 

A- -0.034 

C,D 

Surface  gravity  log  g*  (cgs) 

4.236  ±  0.011 

C,D 

Luminosity  L*  (L©) 

1  QQ+0.24 

0.21 

C,D 

Age  (Gyr) 

3.8  ±  1.0 

C,D 

Planetary  parameters 

Mass  Mp  (Mj ) 

0  669+0  025 
u.ooy_0  030 

A,  B,C,  D 

Radius  Rp  ( Rj ,  equatorial) 

1  99  o+0. 026 

-0.029 

A,  B,C,  D 

Density  pp  (gem-3) 

O  9S9+0.018 
u.  jjz,_0  022 

A,  B,  C,  D 

Surface  gravity  log  gp  (cgs) 

9  074+O.OI6 

7  -0.022 

A,  B,  C 

Orbital  semimajor  axis  a  (AU) 

0  04567+0  00055 

u.uhoo/_0.00046 

E 

Equilibrium  temperature  req  (K) 

1500  ±  200 

F 

Notes. 

A:  based  on  the  photometry. 

B:  based  on  the  radial  velocities  assuming  an  elliptical  orbit. 

C:  based  on  an  SME  analysis  of  the  HIRES  spectra. 

D:  based  on  the  Yale-Yonsei  stellar  evolution  tracks. 

E:  based  on  Newton’s  version  of  Kepler’s  third  law  and  total  mass. 

F:  assumes  Bond  albedo  =  0.1  and  complete  redistribution.  The  uncertainty 
reflects  the  uncertainties  in  these  assumptions. 

detection  of  centroid  shifts  is  on  the  order  of  0. 1  millipixel  or 
somewhat  less.  SeeBatalhaetal.  (2010)  for  a  detailed  discussion 
of  this  approach  for  candidate  vetting. 

3.3.  Precise  Radial  Velocities 

The  precise  velocities  shown  in  Table  1  were  obtained  by 
using  HIRES  on  the  Keck  1  telescope  with  iodine  cell  radial 
velocity  reference.  A  modified  version  of  the  standard  HIRES 
iodine  cell  reduction  pipeline  that  includes  improved  cosmic  ray 
rejection  and  sky  brightness  suppression  was  applied  to  these 
spectra  by  J.  Johnson.  The  reduction  pipeline  modification  was 
necessary  because  of  the  faintness  of  this  star  compared  to  the 
bright  stars  normally  observed  for  radial  velocity  planet  search 
work.  The  nearby  companion  mentioned  earlier  is  sufficiently 
far  away  on  the  sky  that  it  does  not  interfere  with  the  HIRES 
work.  The  phased  velocities  are  shown  in  Figure  2  with  a  fit  to  a 
circular  orbit  whose  phase  and  period  are  fixed  by  the  observed 
transit  times. 

The  Monte  Carlo  analysis  described  below  shows  no  evidence 
of  orbital  eccentricity.  In  this  analysis,  the  eccentricity  and 


Phase 


Figure  2.  Orbital  solution  for  Kepler-6.  The  observed  radial  velocities  obtained 
with  HIRES  on  the  Keck  1  Telescope  are  plotted  in  the  top  panel  together  with 
the  velocity  curve  for  a  circular  orbit  with  the  period  and  time  of  transit  fixed 
by  the  photometric  ephemeris.  The  middle  panel  shows  the  residuals  to  the  fit. 
The  point  potentially  affected  by  the  Rossiter-McLaughlin  effect  is  the  one  at 
the  smallest  phase.  With  this  point  omitted  from  the  fit  the  rms  of  the  residuals 
is  5.7  ms-1.  The  lower  panel  shows  the  line  bisector  variation.  No  structure 
connected  to  the  radial  velocity  curve  is  evident. 

longitude  of  periastron  were  parameterized  as  e  sin  o>  and 
e  cos  co.  The  probability  distribution  for  these  parameters  is 
symmetrical  with  zero  mean  and  a  standard  deviation  of  0.029. 

The  rms  of  the  velocity  residuals  is  5.7ms”1  with  the  point 
nearest  transit  omitted  from  the  fit.  This  point  is  potentially 
affected  by  the  Rossiter-McLaughlin  (R-M)  effect.  Its  residual 
of  12ms”1  is  larger  than  all  the  other  velocity  points  and  is  in 
the  sense  of  a  prograde  orbit.  The  rotational  velocity  given  in 
Table  2  would  result  in  an  R-M  amplitude  of  20  m  s”1 ,  certainly 
detectable  with  further  careful  spectroscopy.  An  analysis  of 
the  line  bisectors  (Figure  2)  shows  an  rms  of  6.8  ms”1  with 
no  organized  structure  that  might  result  from  a  triple  system 
(Mandushev  et  al.  2005). 

4.  DISCUSSION 

A  Monte  Carlo  bootstrap  procedure  developed  by  Rowe  and 
Brown,  outlined  in  Koch  et  al.  (2010b)  and  Borucki  et  al. 
(2010b),  that  simultaneously  fits  the  light  curve  and  the  radial 
velocities,  and  that  is  consistent  with  stellar  evolution  models, 
was  carried  out  to  produce  most  of  the  stellar  and  planetary 
properties  presented  in  Table  2.  The  model  light  curve  (Mandel 
&  Agol  2002)  was  integrated  over  the  long  cadence  integration 
period  prior  to  fitting  to  the  Kepler  data.  We  note  that  the 
nearby  companion  star  alluded  to  earlier  is  only  about  one 
Kepler  pixel  away  from  Kepler-6,  so  its  light  is  included  in 
the  Kepler  aperture.  The  fainter,  more  distant  star  is  partially 
included  in  the  aperture  as  well.  This  has  the  effect  of  diluting 
the  depth  of  the  transit  by  3.3%  ±  0.4%.  This  effect  has  been 
taken  into  account  in  the  modeling  process  including  the  dilution 
uncertainty  due  to  the  partial  inclusion  of  the  more  distant  star. 
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The  derived  planetary  density,  0.352^°0 022  gem-3,  places  it 
in  a  well-populated  area  of  the  mass-radius  relationship  for 
currently  known  extrasolar  planets. 

Several  spectra  obtained  with  HIRES  without  the  iodine  cell 
were  subjected  to  an  SME  analysis  (Valenti  &  Piskunov  1996) 
by  D.  Fischer  to  provide  the  remaining  stellar  properties  reported 
in  Table  2.  In  this  analysis,  log  g  was  frozen  at  the  transit- 
derived  value  since  it  is  more  reliable  than  the  SME  value, 
which  was  4.59  if  allowed  to  float.  The  discrepancy  between 
the  transit  value  of  log  g  and  the  SME  value  is  significant. 
Gilliland  et  al.  (2010b)  and  Nutzman  et  al.  (2010)  find  a  similar 
effect  in  their  comparison  of  log  g  derived  from  transit  fits, 
asteroseismology  constraints,  and  SME  spectral  analysis  in  the 
caseofHD17156  ([Fe/H]  =  +0.24).  In  that  case  the  discrepancy 
in  log  g  is  smaller,  0.12,  but  is  in  the  same  sense.  Bakos  et  al. 
(2009)  find  a  similar  situation  for  HAT-P-13  ([Fe/H]  =  +0.41) 
with  a  log  g  discrepancy  of  +0.14.  This  suggests  that  there 
is  a  systematic  error  in  the  logg  values  derived  from  SME 
analysis  for  high-metallicity  stars.  However,  Burke  et  al.  (2007) 
find  only  a  small  log  g  discrepancy  of  +0.02  in  the  case  of 
XO-2  ([Fe/H]  =  +0.45)  and  the  discrepancy  for  HD80606 
([Fe/H]  =  +0.4)  is  only  +.01  (Naef  et  al.  2001;  Winn  et  al. 
2009).  Further  comparisons  of  spectroscopic  log  g  values  with 
precise  values  derived  from  transits  and  asteroseismology  are 
likely  to  be  fruitful. 
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